from pylab import *

vecs = loadtxt('JLA_WMAP9_SP_vecs.txt')
#vecs = loadtxt('JLA_WMAP9_vecs.txt')

zmax = 1.5
amin = 1./(1.+zmax)
a = linspace(1,amin,20)
z = 1./a-1.

for i in range(20):
    subplot(4,5,i+1)
    plot(z,vecs[i,:],'-.')
#    legend(loc='best',frameon=False)

show()
